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COAGULATION'  IN  TURBULENT  FLOW  - THEORY  AND  EXPERIMENT 


► 


Michael  A.  Uelichatsios^  and  Ronald  F.  Probstein 

Department  of  Mechanical  Engineering,  Massachusetts  Institute 
of  Technology,  Cambridge,  Massachusetts  02139 

Abstract 

The  coagulation  of  colloidal  particles  in  turbulent  flows  is  in- 
vestigated theoretically  and  experimentally.  A coagulation  model  is 
developed  for  destabilized  particles  in  an  isotropic  turbulent  flow. 
Simple  binary  collision  mean  free  path  concepts  are  employed  and  for 
monodisperse  systems  coagulation  rate  relations  are  derived  for  particle 
sizes  less  than  and  larger  than  the  Kolmogorov  microscale  of  turbulence. 
Polydisperse  systems  are  also  considered  and  some  general  results  con- 
cerning their  behavior  are  obtained.  Also  discussed  are  the  effects  of 
interparticle  repulsion  on  the  coagulation  rates. 

Coagulation  experiments  are  reported  on  inside  fully  developed 
turbulent  pipe  flows  using  an  approximately  monodisperse  UCAR  latex 
dispersion  in  which  the  particle  sizes  are  less  than  the  Kolmogorov 
microscale.  The  flow  rate,  destabilizer  concentration  and  volume  frac- 
tion of  the  dispersed  phase  were  varied  in  these  experiments.  The 
experimental  results  for  destabilized  particles  are  shown  to  agree  very 
well  with  the  theoretical  predictions.  Brownian  motion  coagulation 
experiments  for  partially  destabilized  systems  are  compared  with  the 
corresponding  turbulent  motion  experiments  and  the  results  indicate 
that  the  coagulation  efficiency  does  not  appear  to  depend  on  the  par- 
ticle transport  mode. 
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1.  Present  address:  Factory  Mutual  Research,  Norwood,  Massachusetts. 
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NOMENCLATURE 

dimensionless  coefficient  characterizing  polydispersity , Eq.  [15] 
turbulent  friction  coefficient,  Eq.  [29]  (dimensionless) 
particle  diameter  (cm) 
pipe  diameter  (cm) 

particle  volume  probability  density  (cm  ^) 

Eulerian  macroscale  of  turbulence  (cm) 

-3 

particle  number  density  (cm  ) 

-3 

cumulative  particle  number  density  (cm  ) 
flow  rate  in  pipe  (gal/min) 

Reynolds  number,  UD/v  (dimensionless) 
time  (sec) 

turbulent  friction  velocity,  Eq.  [27]  (cm/sec) 

induced  relative  velocity  between  particles  due  to  repulsive 

effects  (cm/sec) 

relative  particle  velocity  during  collision  (cm/sec) 
effective  collision  velocity  for  partially  destabilized 
particles  (cm/sec) 
average  pipe  velocity  (cm/sec) 

3 

particle  volume  (cm  ) 

characteristic  potential  energy  of  interaction  between  particles 
(ergs) 

distance  along  pipe  measured  from  first  sampling  station  (cm) 
Brownian  coagulation  stability  factor  (dimensionless) 
turbulent  coagulation  stability  factor  (dimensionless) 
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Greek 
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vU.t) 
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ui(v,v) 


$(5,5') 


rate  of  turbulent  energy  dissipation  per  unit  mass  (ergs/g  * sec) 
Debye  length  (cm) 

3 1/A 

Kolmogorov  microscale,  (v  /e)  (cm) 

viscosity  of  continuous  phase  (g/cm  * sec) 

2 

kinematic  viscosity  of  continuous  phase  (cm  /sec) 

dimensionless  particle  volume,  n v/ d> 

1/2 

Kolmogorov  time  scale,  (v/e)  (sec) 

volume  fraction  of  dispersed  phase  (dimensionless) 
volumetric  concentration  of  destabilizer  (dimensionless) 
dimensionless  particle  volume  probability  density,  Eq . [13a] 
particle  coagulation  (collision)  rate  (sec  ^) 
retarded  coagulation  rate  (sec 

volume  swept  by  each  particle  per  unit  time  during  collision 

(cm’Vsec) 

dimensionless  collision  frequency  function,  Eq.  [13c] 


Subscripts 

1,2  reference  to  two  monodisperse  systems  1 and  2 with  different 

particle  volume  fractions,  diameters  and  number  densities 
Superscript 


reference  to  particles  of  diameter  d'  and  volume  v' 
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the  coagulation  mechanism  is  of  primary  importance  with  regard  to  the 
behavior,  handling  and  treatment  of  dispersions. 

Particles  dispersed  in  a liquid  medium  generally  acquire  an  elec- 
trical charge,  usually  negative.  As  a result  of  the  potential  field 
surrounding  each  particle,  coagulation  is  inhibited  at  short  inter- 
particle  distances  by  the  repulsive  action  between  the  particles.  It 
follows  that  coagulation  may  be  considered  as  a two-step  process  (1): 
a)  kinetic  transport  of  particles;  b)  actual  particle  aggregation  after 
the  interparticle  repulsion  is  overcome.  Particle  transport  may  take 
place  as  a consequence  of  Brownian  motion,  shear  flow  or  turbulent 
motion.  The  repulsive  interparticle  forces  can  be  neutralized,  that  is, 
the  dispersion  can  be  destabilized,  by  the  addition  of  appropriate 
chemicals  usually  called  coagulants  or  coagulant  aids  (1) . 

It  was  Camp  and  Stein  (2)  who  first  suggested  that  the  coagulation 

rate  in  turbulent  flows  inside  agitator  tanks  can  be  estimated  using  the 

theory  developed  by  Smoluchowski  (3)  for  shear  flow  coagulation.  In 

deriving  their  result  for  the  coagulation  rate  they  replaced  the  rate  of 

strain  in  Smoluchowski 's  expression  by  an  effective  turbulent  rate  of 

1/2 

strain  equal  to  (e/v)  , where  £ is  the  average  turbulent  energy  dissipa- 

tion rate  per  unit  mass  inside  the  agitator  tank  and  v is  the  kinematic 
viscosity  of  the  continuous  phase  [see  also  (1)]. 

Later  Saffman  and  Turner  (4)  proposed  a theory  to  estimate  the 
coagulation  rate  of  drops  in  turbulent  clouds  in  which  they  considered 
the  nature  of  the  flow  field  around  a drop.  The  drop  sizes  in  their 
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theory  were  smaller  than  the  Kolmogorov  microscale  of  turbulence.  The  form 
of  their  results  for  the  collision  rate,  which  they  identified  with  the 
coagulation  rate,  agreed  with  the  equation  proposed  by  Camp  and  Stein  to 
within  a constant  factor. 

Finally  Levich  (5)  developed  a coagulation  model  based  on  a diffusion 
equation  for  the  relative  particle  motion,  and  estimated  the  coagulation 
rate  for  particles  of  a size  less  than  the  Kolmogorov  microscale.  His 
model  led  to  the  same  functional  form  for  the  coagulation  rate  as  the  one 
obtained  by  the  previously  mentioned  authors  but  with  a value  larger  by 
approximately  an  order  of  magnitude.  This  discrepancy  may  be  attributed 
to  the  inability  of  a steady  diffusion  model  to  correctly  describe  the 
relative  motion  of  particles  in  turbulent  flows  due  to  the  non-Markovian 
nature  of  the  turbulent  dispersion  process. 

Experiments  carried  out  in  agitator  tanks  [see  Birkner  and  Morgan 
(6) ] have  verified  the  functional  form  for  the  coagulation  rate  predicted 
by  the  various  theories  but  were  unable  to  establish  which  of  the  pro- 
portionality constants  was  correct  among  the  different  values  predicted. 

The  main  reason  for  the  uncertainty  is  that  the  turbulent  conditions  in 
an  agitator  tank  are  far  from  isotropic  (7),  so  that  the  use  of  average 
turbulent  properties  in  correlating  the  experimental  results  (6)  introduces 
an  error. 

In  this  work  a new  model  is  developed  for  destabilized  systems 
based  on  simple  binary  collision  mean  free  path  concepts  for  noninteract- 
ing particles.  For  monodisperse  systems  coagulation  rate  relations  are 
obtained  for  particle  sizes  less  than  and  larger  than  the  Kolmogorov 
microscale  of  turbulence.  Polydisperse  systems  of  unequal  size  spherical 
particles  are  also  considered  and  some  general  results  concerning  their 
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bchavior  are  obtained.  Also  discussed  are  the  effects  of  interparticle 
repulsion  on  the  rate  of  coagulation  between  particles.  The  theoretical 
results  are  compared  with  coagulation  rate  experiments  carried  out  inside 
fully  developed  turbulent  pipe  flows.  Fully  developed  turbulent  pipe 
flow  was  selected  because  its  characteristics  are  well  known  and  because 
the  flow  is  nearly  isotropic  at  the  core  of  the  pipe  (8) . A further 
reason  for  carrying  out  these  experiments  was  that  no  coagulation  measure- 
ments in  turbulent  pipe  flows  have  been  reported  previously  in  the 
literature. 

KINETIC  MODEL  FOR  COAGULATION  IN  TURBULENT  FLOW 

In  what  follows  the  particles  are  considered  to  be  sufficiently 
small  that  we  may  neglect  inertial  effects  so  that  the  particles  can  be 
assumed  to  follow  the  turbulent  motion,  which  in  our  analysis  we  take  to 
be  isotropic.  Particle-particle  interactions,  gravitational  forces  and 
any  breakup  of  particles  due  to  turbulent  forces  are  neglected.  The 
particle  volume  fraction  is  taken  sufficiently  small  (<t>  'v  3%)  that  only 
binary  collisions  need  be  considered.  At  first  we  examine  a destabilized 
system  in  which  the  interparticle  repulsion  has  been  reduced  to  zero  by 
the  addition  of  appropriate  chemicals  (1) . 

Under  the  assumptions  outlined  we  may  apply  simple  binary  collision 
mean  free  path  concepts  for  noninteracting  particles  to  calculate  the 
coagulation  frequency  associated  with  the  statistical  nature  of  a tur- 
bulent flow.  Thus  for  a raonodisperse  system,  consisting  of  equal  size 
spherical  particles,  the  volume  swept  out  by  each  particle  per  unit  time 
during  collision  will  be 


o!(v, v)  * nd^ 


J 


[1] 
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Here  v is  die  particle  volume,  d its  diameter, Td  the  effective  colli- 
sion cross  section  and  ur  the  relative  particle  velocity  during  collision. 

It  follows  that  the  coagulation  rate,  equal  in  effect  to  the  par- 
ticle collision  rate,  is  given  by 


where  n is  the  particle  number  density  and  the  factor  1/2  ensures  that 
each  collision  of  two  equal  size  particles  is  counted  only  once. 

The  decrease  of  particle  number  density  with  time  due  to  coagula- 
tion will  be 


dn 

dt  * ' 


[31 


The  relative  particle  velocity  appearing  in  Eq.  [1]  is  approximately 
equal  to  the  root  mean  square  relative  turbulent  velocity  between  two 
points  at  a distance  of  a particle  diameter  apart.  From  the  Kolmogorov 
theory  of  isotropic  turbulence  the  relative  velocity  depends  on  the  tur- 
bulent scale  and  its  magnitude  is  given  by  the  relations  (9) 

ur  = (^)1/2  d d < X [4a] 

uf  = 1.37  (ed)^3  d > X [4b] 

ur  * (cL)1/3  d % L [4c] 
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microscale  defined  by  the  relation 


i , 3,  ,1/A 
X - (v  /e) 


where  e is  the  rate  of  turbulent  energy  dissipation  per  unit  mass  and 
v is  the  kinematic  viscosity  of  the  continuous  phase.  We  note  that 
Eq.  [Ab]  is  applicable  only  when  an  inertial  subrange  of  turbulence 
exists,  that  is,  for  sufficiently  large  Reynolds  numbers  where  L >>  X 
(9).  Furthermore,  it  is  assumed  here  that  the  macroscale  is  of  a size 
which  could  be  acquired  by  a dispersed  particle  (e.g.,  L 'v  0.1  cm). 

Using  the  appropriate  value  for  the  relative  velocity  u^  in  Eq . 
[2],  we  get  after  some  rearrangements 


T1  - 0.77 


d < X 


— = A. 11  (h2/3 
$ d 


d > X 


— % 3(r)1/3  4 

$ Ad 


d V L 


In  these  equations  t is  the  Kolmogorov  time  scale  defined  by 


t = (v/e) 


and  (fi  is  the  volume  fraction  of  the  dispersed  phase. 


1 

" 6 n nd 


[8] 
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We  liave  plotted  in  Fig.  1 the  dimensionless  coagulation  rates  ojt/^ 

given  by  Eqs . [6]  as  a function  of  d/A.  In  plotting  Eq . [6c]  we  have 

3 

assumed  by  way  of  example  that  L/A  = 10  , which  implies  sufficiently 
large  Reynolds  numbers  that  there  exists  an  inertial  subrange  of  tur- 
bulence. In  Table  I the  theoretical  dimensionless  coagulation  rate  of 
Eq.  [6a]  derived  for  particle  sizes  less  than  the  Kolmogorov  microscale 
is  compared  with  the  comparable  results  of  previous  authors.  The  value 
of  9.24  attributed  to  Levich  was  obtained  by  taking  the  constant  of 
proportionality  in  his  effective  turbulent  diffusion  coefficient  to  be 
equal  to  /1/15,  consistent  with  the  Kolmogorov  result  (cf.  Eq . [4a1). 

TABLE  I 

THEORETICAL  COAGULATION  RATES  FOR  PARTICLE  SIZES  SMALLER  THAN 
THE  KOLMOGOROV  MICROSCALE 


Author 

(i)T/$ 

Camp  and  Stein  (2) 

1.27 

Saffman  and  Turner 

(M 

1.23 

Levich  (5) 

9.24 

Present  study  (Eq. 

[6a]) 

0.77 

For  a spatially  homogeneous  polydisperse  system,  consisting  of  un- 
equal spherical  particles  of  volume  v,  it  is  necessary  to  introduce  the 
particle  volume  probability  density  f(v,t)  at  a time  t.  The  kinetic 
equation  governing  the  temporal  evolution  of  f or  population  balance  is 
then  written  (10) 
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df  (v,t) 
dt 


w(v',v-v')  f(v\t)  f(v-v',t)  dv' 


- | uJ(v,v')  f(v,t)  f(v',t)  dv'  [9] 

0 

where  w(v,v')  is  the  volume  swept  per  unit  time  by  each  of  the  particles 
with  volumes  v and  v'.  According  to  our  simple  collision  model, 

w(v.v')  - TT(d-  y-d-'-)2  ur  [10] 

where  d,d'  are  the  diameters  of  particles  of  volume  v and  v'  respec- 
tively, and  the  relative  velocity  uf  is  equal  to  the  root  mean  square 
velocity  between  two  points  at  a distance  (d  + d')/2  in  an  isotropic 
turbulent  flow  [see  (9) ] . 

Integrating  Eq.  [9]  with  respect  to  v over  all  volumes  and  introduc- 
ing the  definitions  of  the  particle  number  density  n and  volume  fraction  <t> 


n = j f(v,t)  dv  [11a] 

0 

oo 

0 = j vf(v,t)  dv  [lib] 

0 

leads  to  the  kinetic  equation 

OO  oo 

~ = - wn  [ j | *U',t)  dc  d£ ' ] [12] 

0 0 


in  terms  of  the  dimensionless  variables  (10) 
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1 “ —■ 2 f(v,t) 
n 


[13a] 


U3b] 


&(£,*')  - k f 2S(v,v') 

Z Co 


[13c] 


We  note  that  for  a monodisperse  system  with  the  same  volume  fraction 
w in  Eq.  [12]  is  the  coagulation  rate  defined  in  Eqs.  [2]  and  [6],  if 
we  assume  that  no  more  than  one  value  of  ur  in  Eq . [4]  is  applicable 
for  the  entire  spectrum  of  particle  sizes  in  evaluating  the  volume 
swept  per  unit  time  by  Eq . [10]. 

For  dispersed  systems  where  the  particle  sizes  are  less  than  the 
Kolmogorov  microscale  we  can  readily  show  from  Eqs.  [4a]  and  [10]  that 


2s  «;,£') 


I(,l/3+c,l/3)3 


[14] 


It  follows  that  we  may  write  Eq.  [12]  in  the  form 


dn 

dt 


- wn  A(t) 


[15a] 


with  A(t)  a coefficient  characterizing  the  polydisperseness  of  the  system 
defined  by 


A(t) 


00  oo 


0 0 


S’1/3)3 


MZ.t)  <H5\t)  d£' 


[15b] 


and  with  w given  by  Eq.  [6a]. 
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The  coefficient  A(t)  will  be  time  dependent  since  as  recently  shown 
by  Pulvermacher  and  Ruckenstein  (11)  a self  similar  solution  for  iH£,t) 
independent  of  time  cannot  exist  for  the  collision  frequency  function 
of  Eq.  [14). 

Another  important  result  we  can  derive  is  that  A(t)  satisfies  the 
inequalities 

0.25  S A(t)  < 1 [16] 

The  above  result  may  be  derived  from  the  definition  of  A(t)  in  Eq. 

[15b]  through  the  use  of  the  algebraic  inequalities 

U + 5')  < UU 3 + S,1/3)3  * 4(5  + 5')  [17] 

and  the  conditions  obtained  from  introducing  the  dimensionless  variables 
of  Eq.  [13]  into  the  definitions  of  Eq.  [11], 

00 

| <|»(5,t)  d£  » 1 [18a] 

0 

ao 

1 CiKC.t)  d5  - 1 U8b] 

0 

The  coefficient  A(t)  can  be  shown  to  take  on  the  value  of  1 for  a 
monodisperse  system,  and  the  minimum  value  of  0.25  for  a system  con- 
sisting of  two  monodisperse  systems  with  widely  differing  volume 
fractions  (<t>^  <<  <f>2)  where  the  number  densities  are  also  widely  dif- 
ferent but  in  the  inverse  order  (n^  >>  n2) . The  derivation  of  this 
result  is  outlined  in  the  Appendix. 
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EFFECTS  OF  INTERPARTICLE  REPULSION 

In  the  analysis  so  far  it  has  been  assumed  that  charge  repulsion 
effects  between  particles  could  be  neutralized  by  the  addition  of  appro- 
priate chemicals  (1) . Colloidal  particles  usually  acquire  a negative 
surface  ciiarge  in  solution  so  that  they  tend  to  repel  each  other  upon 
collision.  Thus  coagulation  between  particles  is  retarded  or  totally 
inhibited,  when  the  particle  surface  potential  (zeta  potential)  is 
different  from  zero. 

Due  to  the  complex  character  of  turbulent  flows  it  is  not  feasible 
to  describe  analytically  the  retarding  effects  on  coagulation  of  inter- 
particle potentials.  In  this  regard,  the  use  of  the  turbulent  diffusion 
model  suggested  by  Levich  (5),  in  which  potential  effects  would  be  incor- 
porated as  was  done  by  Fuchs  (12)  for  Brownian  motion,  is  not  justified 
because  of  the  non-Markov ian  nature  of  turbulent  flows. 

In  what  follows  a rather  qualitative  description  of  the  effects  of 
interparticle  repulsion  is  proposed  for  particles  of  a size  less  than 
the  Kolmogorov  microscale.  The  characteristic  potential  energy  of  inter- 
action between  the  particles  is  designated  by  V and  typically  may  be 
taken  equal  to  the  maximum  interparticle  potential  energy.  The  range 
over  which  this  potential  energy  is  effective  is  characterized  by  the 
Debye  length  n which  is  the  characteristic  thickness  of  the  double  layer 

surrounding  the  colloidal  particles  (13).  In  general,  the  particle  size 

o 

is  large  compared  to  the  Debye  length  n,  which  typically  is  like  50  A. 

The  other  force  acting  on  the  particles  is  fluid  dynamic  in  origin 
and,  for  particle  sizes  less  than  the  Kolmogorov  microscale,  is  the  low 
Reynolds  number  Stokes  drag 
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F % 3 irp  u^d  [19] 

where  uf  is  the  relative  particle  velocity  during  collision  (cf.  Eq.  [1]) 
and  u is  the  viscosity  of  the  continuous  phase.  The  approximately  equal 
sign  has  been  used  here  because  the  drag  force  will  be  somewhat  smaller 
than  indicated  as  a result  of  the  fluid  dynamic  interactions  between  the 
particles  as  they  approach  each  other. 

Comparing  the  electrostatic  potential  force  to  the  fluid  dynamic 
drag  force,  it  can  be  seen  from  dimensional  arguments  that  the  retarded 
coagulation  rate  w*  is  expressible  through  the  functional  relation 

in*  ■ w fen  ( — ) [20] 

* u urd  n 

Here  w is  the  coagulation  rate  for  the  fully  destabilized  system  and  the 
function  represents  a coagulation  efficiency  which  tends  to  1 as  V + 0 
and  to  0 as  V + *. 

The  motion  of  a particle  of  a size  less  than  the  Kolmogorov  micro- 
scale in  the  potential  field  of  another  particle  is  characterized  by 

three  time  scales,  the  particle  motion  characteristic  relaxation  time 

2 1/2 

(d  / 18 v) , the  characteristic  Kolmogorov  eddy  time  scale  [(v/e)  ],  and 

the  time  required  for  the  particle  to  move  a distance  equal  to  the 
Debye  length  (n/u^).  We  may  estimate  that  the  relaxation  time  for 
micron  size  particles  is  much  less  than  the  time  required  for  the  par- 
ticle to  move  a Debye  length,  while  the  characteristic  turbulent  eddy 
time  is  much  larger.  It  follows  that  the  change  in  turbulent  velocity 


-15- 


during  the  time  the  particles  interact  can  be  neglected.  The  amount  by 
which  the  coagulation  frequency  is  reduced  when  the  interparticle 
potential  energy  V is  small  may  then  be  estimated  in  the  following  way. 
At  a distance  of  the  order  of  a Debye  length  n,  the  repulsive  potential 
energy  between  the  particles  induces  a relative  particle  velocity  of 
the  order  of 


ui  % 3ttu  d 


[21] 


Therefore  at  short  interparticle  distances  the  instantaneous  relative 
turbulent  particle  velocity  will  be  reduced  by  an  amount  given  by  u^ . 

The  effective  collision  velocity  between  the  particles  can  be 
found  from  the  requirement  that  the  instantaneous  relative  turbulent 
velocity  must  be  absolutely  larger  than  the  induced  velocity  in  order 
for  a collision  to  lead  to  actual  particle  coagulation.  For  small 
values  of  the  induced  velocity  in  comparison  with  the  root  mean  square 
velocity  difference  uf  (cf.  Eq.  [1]),  the  effective  collision  velocity 
may  be  approximated  by 


^r^eff 


u 

r 


[22] 


Using  the  same  simple  mean  free  path  model  defined  through  Eq . [2], 
with  an  interparticle  velocity  (ur)e£f>  wc  obtain  the  following 
expression  for  the  reduced  coagulation  rate 
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wliere  we  have  used  the  value  of  u^  given  by  Eq . [21]  and  have  assumed 

n/d  <<  1. 

EXPERIMENTS 

Experiments  were  carried  out  to  measure  the  coagulation  rate  of 
colloidal  particles  inside  a fully  developed  turbulent  pipe  flow.  The 
principal  advantage  of  this  configuration  in  contrast  to  say  a stirred 
tank,  is  that  the  turbulent  characteristics  of  such  a flow  are  well 
known  [see,  e.g.,  (8)].  Furthermore,  in  a fully  developed  flow  the 
turbulence  is  nearly  isotropic  in  the  pipe  core.  Commercial  latex  par- 
ticles having  a narrow  size  distribution  were  dispersed  in  water.  The 
particle  sizes  and  the  number  of  particles  per  unit  volume  were  measured 
using  a Coulter  counter  [see,  e.g.,  (14)]. 

Figure  2 shows  a diagram  of  the  experimental  setup.  The  flow  is 
established  by  means  of  a centrifugal  pump.  The  flow  rate  is  measured 
by  a rotameter  and  controlled  by  the  overflow  arrangement  shown  in  the 
figure.  Two  honeycomb  filters  are  used  in  the  feed  line  to  filter  the 
supply  water.  The  stable  latex  emulsion  was  released,  at  a measured 
rate,  at  station  I into  the  pipe  flow  from  a pressurized  tank,  through 
a solenoid  valve  controlled  by  an  electric  timer.  Destabilizer  is  in- 
jected, at  a measured  rate,  at  station  II,  where  the  latex  particles 
are  known  to  be  completely  mixed  with  the  turbulent  flow.  This  in- 
jection system  is  similar  to  that  for  the  emulsion. 

Latex  particles  manufactured  by  Union  Carbide  (878  UCAR)  with  an 
average  size  of  0.6  um  were  used  in  these  experiments.  A mixture  of 
1 mol/liter  of  HC1  and  5.6  mo] /liter  of  CaC^  was  used  as  the  de- 


stabilizer. A special  construction  for  the  injection  port  at  station 
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II  ensured  complete  mixing  of  the  destabilizer  with  the  flow  in  a short 
distance.  Sufficient  pipe  length  is  provided  after  each  bend,  so  that 
fully  developed  turbulent  flow  prevails. 

After  steady  flow  conditions  are  established  (for  the  present 
experiments,  approximately  10  sec  after  the  release  of  the  latex  emul- 
sion) samples  are  withdrawn  at  five  stations  along  the  pipe.  The  sample 
at  each  station  is  withdrawn  automatically  by  a solenoid  valve,  through 
a 1/4  inch  I.D.  pipette  in  the  wall  of  the  pipe.  The  velocity  through 
the  pipette  was  much  larger  ('v  500  cm/sec)  than  the  average  pipe 
velocity  (-v  180  cm/sec),  so  that  the  withdrawn  sample  is  representative 
of  the  conditions  at  the  core  of  the  pipe  flow.  The  coagulation  inside 
this  pipette  was  estimated  to  amount  to  less  than  27.  of  the  total 
particle  number  density.  As  shown  in  Fig.  1 each  sample  is  diluted 
directly  in  a 1 liter  flask  containing  distilled  water.  The  dilution 
is  sufficient  to  ensure  that  Brownian  coagulation  effects  are  negligible, 
prior  to  the  time  at  which  the  actual  analysis  of  each  sample  is  carried 
out . 

Prior  to  carrying  out  the  measurements,  each  sample  is  further 
diluted  by  a known  amount  in  a pure  aqueous  solution  of  1%  sodium 
chloride.  These  samples  are  then  analyzed  using  the  Coulter  counter. 

The  theory  behind  the  operation  of  the  Coulter  counter  is  quite  simple. 
The  dispersed  particles  suspended  in  the  1%  aqueous  solution  of  sodium 
chloride  are  forced  to  flow  through  a small  aperture  which  has  an  im- 
mersed electrode  on  each  side.  As  each  particle  passes  through  the 
aperture  it  replaces  a volume  of  electrolyte  equal  to  its  own  value 
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within  the  aperture,  thus  momentarily  changing  the  resistance  between 
the  electrodes.  This  produces  a voltage  pulse  of  short  duration  which 
has  a magnitude  proportional  to  the  particle  volume.  The  resulting 
series  of  pulses  is  electronically  amplified,  rated  and  counted.  Thus 
for  each  sample  not  only  is  the  total  number  of  particles  in  a measured 
volume  counted,  but  the  particle  size  distribution  is  also  measured. 

Apertures  of  different  sizes  can  be  used  for  different  ranges  of  par- 
ticle size.  Detailed  information  on  this  measurement  technique  may  be 
found  in  the  previously  cited  paper  of  Wachtel  and  La  Mer  (14). 

Calibration  of  the  counter  was  done  using  latex  particles  of  a 
known  size.  Sufficiently  diluted  samples  (<t  ^ 10  were  used  for  the 
analysis  so  that  corrections  due  to  coincidental  passing  of  two  or  more 
particles  through  the  aperture  were  unnecessary.  Aperture  sizes  of  19 
um,  50  pm  and  200  pm  were  used  in  the  present  experiments. 

Additional  details  of  the  apparatus,  including  dimensions  of  the 
various  components,  may  be  found  in  (15). 

RESULTS 

Prior  to  carrying  out  detailed  measurements  in  the  turbulent  pipe 

flow  system  described,  the  Brownian  coagulation  rate  was  measured  for 

different  destabilizer  concentrations  for  the  latex  dispersion  to  be 

used  in  the  pipe  flow  experiments.  Samples  from  a dilute  dispersion 
-3  -4 

(it  ■ 10  to  10  ) were  taken  from  a beaker  at  specified  time  intervals 

and  analyzed  with  the  Coulter  counter.  The  purpose  of  these  tests  was 
to  determine  the  destabilizer  concentration  required  for  the  dispersion 
to  be  fully  destabilized. 

I 
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Figure  3 shows  the  experimental  results.  In  this  figure  the 
stability  factor  is  plotted  as  a function  of  the  volumetric  de- 
stabilizer concentration  The  stability  factor  is  the  ratio  of  the 

coagulation  rate  for  a fully  destabilized  system  to  the  coagulation  rate 
for  a partially  destabilized  system  (13)  (cf.  Eq . [20]).  From  this 
figure  we  conclude  that  the  destabilizer  concentration  necessary  for 

'■ 

complete  destabilization  of  the  present  dispersion  is  approximately  A. 5%. 
The  experimental  linear  fit  for  destabilizer  concentrations  less  than 
4.5%  has  been  drawn  consistent  with  previous  results  (16). 

Experiments  in  the  turbulent  flow  system  shown  in  Fig.  2 were 
carried  out  in  which  the  flow  rate  through  the  pipe,  the  volume  fraction 

of  the  dispersed  phase  and  the  destabilizer  volume  concentration  were 

t 

varied.  The  flow  rates  Q ranged  between  5 and  15  gal/min  for  which  the 

corresponding  Reynolds  number  based  on  pipe  diameter  ran  between  17,000 

and  51,000.  The  volume  fraction  $ of  the  latex  particles  was  varied 
-3  -3 

between  10  and  6 x 10  while  the  volumetric  fraction  of  the  de- 
stabilizer was  controlled  between  1 and  10%.  Brownian  coagulation 
was  always  negligible  compared  to  turbulent  coagulation  and  the  particle 
size  in  these  experiments  was  always  less  than  the  Kolmogorov  microscale. 

In  each  experiment  the  particle  size  distribution  was  measured  at 
every  sampling  station  along  the  pipe  (see  Fig.  2) . A typical  set  of 
results  is  shown  in  Fig.  4 for  a pipe  flow  rate  of  15  gal/min.  In  this, 
figure  the  cumulative  particle  number  density  n^  is  plotted  as  a func- 
tion of  the  particle  volume  at  different  sampling  stations.  The 
cumulative  particle  number  density  gives  the  number  of  particles  per 
unit  volume  of  the  continuous  phase  with  volumes  larger  than  a certain 
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volume  v.  It  is  defined  by  (cf.  Eq.  [11a]) 


n 

v 


j f(v, t)  dv 
v 


[24] 


From  Fig.  4 it  can  be  seen  that  the  total  particle  number  density 
decreases  along  the  pipe  due  to  coagulation. 

Figure  5 shows  the  same  experimental  points  as  in  Fig.  4 plotted 
in  the  reduced  dimensionless  similarity  variables  suggested  by  Eqs . 
[13].  This  figure  would  seem  to  Indicate  that  over  the  length  of  pipe 
in  which  the  experiments  were  conducted  the  cumulative  particle  size 
distribution  at  each  sampling  station  remains  approximately  similar  to 
the  initial  particle  size  distribution  at  the  first  station.  A best 
fit  curve  through  these  points  is  drawn  in  Fig.  5.  Using  this  curve 
the  coefficient  A(t) , which  was  introduced  in  Eq.  [15]  to  characterize 
the  polydisperseness  of  the  system,  is  found  to  be 


A(t)  - 0.90 


[25] 


where  any  time  dependence  has  been  suppressed  through  the  use  of  the 
"best  fit  similarity  curve." 

Due  to  the  limited  length  of  the  pipe  in  the  present  setup,  no 
experiments  could  be  run  for  times  much  larger  than  the  characteristic 
coagulation  time,  in  order  to  test  whether  the  particle  size  distribu- 
tion was  or  was  not  self  similar.  However,  as  we  have  discussed 
previously,  a self  similar  solution  independent  of  time  cannot  exist 
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for  tl>e  collision  frequency  function  consistent  witli  the  present 
theoretical  model  (Eq.  [14]). 

Figure  6 is  a typical  log-linear  plot  of  the  total  particle  number 
density  as  a function  of  distance  along  the  pipe.  The  ordinate  in  this 
figure  is  the  ratio  of  the  total  particle  number  density  at  a sampling 
station  to  the  total  particle  number  density  at  station  one.  The 
abscissa  is  the  time  required  for  the  dispersion  to  flow  to  the  sampl- 
ing station 

t =•  x/U  [26] 

where  x is  the  distance  along  the  pipe  measured  from  the  first  sampl- 
ing station  and  U is  the  average  pipe  velocity. 

For  the  same  degree  of  destabilization,  the  experimental  points 
lie  on  a straight  line  in  agreement  with  the  theoretical  Eq . [3]  for 
a constant  coagulation  rate  w.  As  may  be  seen  in  the  figure,  the  slope 
of  the  lines  which  define  the  coagulation  rate  increases  absolutely 
with  an  increase  in  the  volumetric  destabilizer  concentration 

In  Fig.  7 we  have  shown  a typical  set  of  results  for  the  measured 
coagulation  rates  as  a function  of  destabilizer  concentration  V The 
Kolmogorov  time  scale  (Eq.  [7])  was  calculated  from  the  characteristics 
of  turbulent  pipe  flow.  Experiments  by  Laufer  (8)  have  shown  that 
fully  developed  turbulent  pipe  flow  is  very  close  to  homogeneous  and 
isotropic  to  within  about  90%  of  the  pipe  radius.  As  previously 
described,  the  sampling  stations  in  the  experiment  were  designed  in  such 
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a way  t.iat  the  sample  withdrawn  from  the  pipe  is  representative  of  the 
conditions  at  the  core  of  the  pipe  flow.  It  follows  that  in  calculat- 
ing the  Kolmogorov  time  scale,  the  turbulent  characteristics  repre- 
sentative of  the  core  of  the  pipe  flow  are  the  ones  to  be  used. 

From  the  experimental  work  of  Laufer  (3)  the  rate  of  energy  dissipa- 
tion per  unit  mass  corresponding  to  the  core  of  the  pipe  flow  is 


3 


Here  L)  is  the  pipe  diameter  and  u^  the  pipe  friction  velocity,  which 
may  be  found  from  the  relation  (17) 


u^  ■ U /c j/ 8 


[28] 


where  c^  is  the  turbulent  friction  coefficient.  Measurements  of  the 
pressure  drop  in  the  piping  system  of  the  experiments  gave  values  of 
Cj  equal  to  those  obtained  from  the  standard  friction  coefficient 
for  turbulent  flow  in  smooth  pipes  (17), 


Cf  “ 


0.316 


Re 


1/A 


[29] 


where  Re  is  the  pipe  flow  Reynolds  number  (Re  = UD/v). 

It  can  be  seen  from  Fig.  7 that  for  <t>D  ^ 5%  the  reduced  coagula- 
tion rate  wt/<| > % 0.78  and  thereafter  remains  essentially  constant 
with  increasing  destabilizer  concentration.  This  indicates  that  the 
dispersion  is  effectively  fully  destabilized  for  destabilizer 


concentrations  larger  than  5%.  Taking  into  account  the  correction  factor 
characterizing  the  polydispers ity  of  the  system  (Eq.  [25])  the  reduced 
coagulation  rate  for  a monodisperse  system  becomes  (see  Eq.  [15a]), 


o)T  _ 0.78 
0 “ 0.90 


0.86 


[30] 


It  can  be  seen  that  this  value  agrees  very  well  with  the  theoretical 
result  of  Eq . [ 6a] . 

For  destabilizer  concentrations  less  than  5%  a linear  fit  has  been 
drawn  through  the  experimental  points  of  Fig.  7,  in  analogy  with  the 
behavior  observed  in  Brownian  coagulation.  However,  we  note  that  this 
behavior  is  in  agreement  with  the  predicted  dependence  shown  in  Eq. 

[23] . This  follows  from  the  fact  that  near  the  critical  destabilizer 
concentration  for  full  destabilization  of  the  dispersion  the  potential 
energy  V is  linear  in  In  [see  (18)]. 

In  Fig.  8 we  have  plotted  the  turbulent  stability  factor  W 
as  a function  of  destabilizer  concentration  for  different  flow  rates. 

The  stability  factor  is  the  ratio  of  the  coagulation  rate  for  the  fully 
destabilized  system  to  the  coagulation  rate  for  the  partially  destabilized 
system  and  is  just  the  inverse  of  the  "coagulation  efficiency"  function 
defined  in  Eq.  [20],  We  may  make  the  following  remarks  concerning  the 
results  of  Fig.  8: 

a.  The  slope  of  the  linear  fit  for  partially  destabilized 
systems  (<$>^  ^ 5%)  is  independent  of  the  flow  rate  and  hence  of  the 
turbulent  rate  of  energy  dissipation  per  unit  mass  (Eq.  [27]).  This 
does  not  agree  with  the  approximate  result  represented  by  Eq.  [23], 
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where  the  relative  velocity  is  a function  of  the  energy  dissipation 
c (see  Eq . [4a]) . 

b.  The  slope  of  the  linear  fit  in  Fig.  3 is  almost  the  same 
as  the  slope  for  the  Brownian  stability  factor  in  Fig.  3,  although  we 
might  expect  that  turbulent  motion  would  lead  to  higher  relative  rates 
of  coagulation  when  compared  with  Brownian  motion.  This  same  behavior 
has  been  observed  by  other  investigators  in  experiments  conducted  in 
agitator  tanks  [see  (19)  J. 

c.  : « --er.ta.  value  of  the  destabilizer  concentration 

for  full  aesta:  : ..  :.  ;-ersion  is  approximately  the  same 

for  Brownian  a:.,  turbulent  coagulation  (5%). 

In  the  course  of  tne  present  study  experiments  were  carried  out 
with  particles  larger  than  the  Kolmogorov  microscale.  For  this  purpose 
the  UCAR  latex  emulsion  was  allowed  to  coagulate  in  the  emulsion  tank 
(see  Fig.  2)  to  the  required  degree  by  the  addition  of  small  amounts 
of  destabilizer.  The  same  procedure,  as  described  in  the  previous 
section,  was  used  to  carry  out  an  experimental  run.  Unfortunately,  the 
results  could  not  be  intelligently  interpreted  because  breakup  of  the 
particles  took  place  simultaneously  with  coagulation.  A comprehensive 
study  and  understanding  of  the  breakup  rate  of  particles  in  turbulent 
flows  is  not  available  at  this  time  to  enable  the  proper  interpretation 
of  the  results  of  turbulent  pipe  flow  coagulation  experiments  with 
particles  larger  than  the  Kolmogorov  microscale. 

CONCLUSIONS 

A model  based  on  simple  binary  collision  mean  free  path  concepts 
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was  developed  for  the  prediction  of  the  coagulation  rates  in  isotropic 
turbulent  flows  of  destabilized  monodisperse  particles  less  than  and 
larger  than  the  Kolmogorov  microscale.  The  same  simple  mean  free  path 
concepts  were  applied  to  the  description  of  the  temporal  evolution  of 
the  particle  volume  probability  density  in  coagulating  polydisperse 
systems.  It  was  shown  that  polydisperseness  of  the  system  resulted  in 
a decrease  in  the  total  particle  number  density  coagulation  rate  in 
comparison  with  the  rate  for  a monodisperse  system  of  the  same  volume 
fraction.  For  particles  of  a size  less  than  the  Kolmogorov  microscale 
an  estimate  was  made  of  the  reduction  in  the  coagulation  rate  as  a 
result  of  the  interparticle  repulsion  associated  with  charged  particles 
in  a partially  destabilized  system.  This  estimate  was  based  on  com- 
parison of  the  fluid  mechanical  (Stokes  drag)  forces  to  the  electro- 
static potential  forces. 

Experimental  results  on  coagulation  of  destabilized  particles  in 
fully  developed  turbulent  pipe  flows,  where  the  particle  sizes  are 
small  compared  to  the  Kolmogorov  microscale,  showed  excellent  agree- 
ment with  the  theoretical  prediction  of  the  coagulation  rate.  The 
coagulation  efficiency  for  partially  destabilized  systems  appeared  to 
be  independent  of  the  particle  transport  mode,  that  is,  Brownian  or 
turbulent,  as  has  been  observed  by  previous  investigators  (19) . 
Experiments  in  which  the  dispersed  particle  sizes  were  larger  than  the 
Kolmogorov  microscale  could  not  be  successfully  analyzed  because  of 
simultaneous  breakup  of  the  particles  in  the  flow. 
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APPENDIX 

In  a oionodisperse  system  the  dimensionless  particle  size  distribu- 
tion defined  by  Eq.  [13a]  Is  a delta  function  centered  at  1.  It  follows 
directly  from  Eq.  [15b]  that  for  a monodisperse  system  at  the  onset  of 
coagulation  the  coefficient  A(t)  - 1. 

In  a system  consisting  of  two  monodisperse  systems  with  particle 
volume  fractions,  diameters  and  number  densities  d^,  n^  and  <p^,  d^, 
n^  respectively,  the  evolution  of  the  total  particle  density  with  time 
at  the  onset  of  coagulation  is  given  by 


d^  + n2) 

dt 


“ 2 “ll  nl 


“12  nln2  2 “22n2 


[A.l] 


Here  the  symmetric  tensor  ■ 1,2)  represents  the  volume  swept 

per  unit  time  by  particles  with  diameters  and  d^  during  collision 
and  can  be  found  from  Eq.  [10].  For  particle  sizes  less  than  the 
Kolmogorov  microscale 


,di  + dj,3 
“ij  " °( 2' ' 


[A. 2] 


where  the  proportionality  coefficient 


/I  ,e.  1/2 

a - * / -rr  (~) 


15  'vy 


[A. 3] 


Then  Eq . [A.l]  can  be  rewritten  in  the  following  form 
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d(”1  + n2) 


(n^  + n2> 


4(njd3  + n2d2)  + (dx  + d2)3  n^ 
4(n^  + n2>  (n^d^  + n2d2) 


[A.4] 


wherein  Eqs.  [A. 2]  and  [A. 3]  have  been  used  and  w is  the  coagulation 
rate  of  a monodisperse  system  of  volume  fraction  (<t>^  + <f>2)  (cf.  Eq . 


[15a]). 


It  follows  that  the  coefficient  A(t=0)  (cf.  Eq.  [15a]  at  the  onset 


of  coagulation)  will  be  equal  to 


A(t-O) 


4(n^d3  + n2d3)  + (dx  + d2)3  n^ 
4^  + n2)  (n^d3  + n2d3) 


[ A.  5 ] 


After  some  rearrangements  it  is  found  that 


...  ...  . 3(dl-d2>  <dl+d2>nln2 

A(t-O)  -1-T 3 j— 

(nx  + n2)  (n1d1  + n2d2) 


[A. 6] 


For  the  special  case  of  two  widely  differing  volume  fractions  (<J>^  <<  4>2) 
and  widely  differing  number  densities  in  the  inverse  order  (n^  >>  n7) 
we  have  using  Eq.  [8] 


nl  d2  3 

1 «zr  « (/> 

n2  dl 


[A. 7] 


and  the  value  of  the  coefficient  A(t)  becomes 


A(t-O)  - 0.25 


[A.  8] 
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Theoretical  reduced  coagulation  rate  between  particles  in  isotropic  turbulent  flows  as 
function  of  the  ratio  of  the  particle  size  (d)  to  the  Kolmogorov  microscale  (X). 
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Fig.  6.  Particle  number  density  relative  to  value  at  first  sampling 
station  as  a function  of  distance  along  pipe  measured  in 
time  for  dispersion  to  flow  to  sampling  station. 


Fig.  8.  Stability  factor  as  a function  of  destabilizer  concentration 

for  turbulent  coagulation. 


